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Near a bifurcation point a system experiences critical slowing down. This leads to scaling behavior 
of fluctuations. We find that a periodically driven system may display three scaling regimes and 
scaling crossovers near a saddle-node bifurcation where a metastable state disappears. The rate of 
activated escape W scales with the driving field amplitude A as \nW {Ac — ^4)*, where Ac is 
the bifurcational value of A. With increasing field frequency the critical exponent ^ changes from 
^ = 3/2 for stationary systems to a dynamical value ^ = 2 and then again to ^ = 3/2. The analytical 
results are in agreement with the results of asymptotic calculations in the scaling region. Numerical 
calculations and simulations for a model system support the theory. 
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I. INTRODUCTION 

Thermally activated transitions are at root of many 
physical phenomena: diffusion in solids, protein folding, 
and nucleation are examples. It is important to under- 
stand how transitions occur, particularly in systems away 
from thermal equilibrium. Full understanding would in- 
clude a description of the underlying dynamics and the 
transition probabilities. Owing to their exponential sen- 
sitivity, these probabilities provide an important means 
of characterizing a system. However, in many cases acti- 
vation barriers are high which leads to very low transition 
rates and impedes precise experimental studies. 

The barrier for escape from a metastable state is re- 
duced when the system is close to a bifurcation (critical, 
or spinodal) point where the state disappears. For sys- 
tems that display hysteresis such a bifurcation point cor- 
responds to the switching point on the hysteresis loop. 
The idea of bringing the system close to the bifurcation 
point (l| has been used in studying activated switching 
in Josephson junctions where it has become a 

standard technique for determining the critical current. 
This idea is also used in studies of activated magnetiza- 
tion reversals in nanomagnets 0, 0, II] ■ 

Experiments on nanomagnets and Josephson junctions 
are often performed by ramping the control parameter 
(magnetic field or current) and measuring time distribu- 
tion of escape events 0. In interpreting the data it is 
usually assumed that, for sufficiently slow ramp rates, 
the system remains quasistationary. In this approxima- 
tion the barrier height, i.e., the activation energy of a 
transition i?, usually scales with the control parameter 
?7, measured from its critical (bifurcational) value rjc = 0, 

as 773/2 

ScaUng of R near a bifurcation point is related to slow- 
ing down of one of the motions (llj, i.e., the onset of a 
"soft mode" . The relaxation time of the system di- 
verges as the control parameter — > 0. Therefore if 
depends on time, even where this dependence is slow the 
assumption of quasistationary may become inapplicable 
for small 77. 

In this paper a theory of activated transitions is devel- 



oped for periodically modulated systems. In such sys- 
tems the notion of a stable state is well-defined irre- 
spective of the modulation rate, and the applicability of 
the quasistationary approximation can be carefully stud- 
ied. It turns out that, unexpectedly, near a critical point 
this approximation breaks down even where the relax- 
ation time t-c is still much smaller than the driving period 

We show that an interplay between the critical slow- 
ing down and the slowness of time-dependent modula- 
tion leads to a rich scaling behavior of the transition rate 
and to crossovers between different scaling regions. This 
behavior near a bifurcation point is system-independent 
and has no counterparts in stationary systems. We find 
three regions in which the activation energy scales as 
R (X ij^. As the parameters change, for example with 
the increase of the modulation frequency lof, the critical 
exponent ^ varies from 3/2 to 2 and then again to 3/2. 
Our numerical calculations and Monte Carlo simulations 
for a model system agree with the general results. 

Activated transitions in periodically driven systems 
were investigated earlier in various contexts |l2l ll^ 
II1II1.17, 18, 19, 20], stochastic resonance and diffusion 
in modulated ratchets being recent examples |2lLl2^l23| . 
In this paper we study the previously unexplored region 
of driving amplitudes close to critical and reveal the uni- 
versality that emerges. 

A qualitative picture of motion near a bifurcation point 
can be obtained if one thinks of the system as a particle 
in a potential U{q,t) that oscillates in time with period 
Tp, see Fig. ^ Such particle has periodic stable and 
unstable states, qa{t) and qb{t)- In the adiabatic limit 
LUp —I- they lie at the minimum and local maximum of 
the potential in Fig. As the modulation amplitude 
A increases, the states become close to each other for a 
portion of the period Tp, see Fig.^. The barrier height 
reaches its minimum during this time, and this is when 
the system is most likely to escape from the potential 
well. 

As A further increases, the states qa,b{t) become dis- 
torted to avoid crossing, see Fig.^, and the adiabatic ap- 
proximation becomes inapplicable. The parameter range 
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FIG. 1: (a) An oscillating potential barrier. In the limit of 
slow modulation, the stable and unstable periodic states qa 
and qt are the instantaneous positions of the potential mini- 
mum and barrier top, respectively, (b) For slow modulation, 
when the driving amplitude A is close to A^"^, the states qa,b{t) 
come close to each other once per period, (c) As A further 
approaches Ac, the states qa.bit) become skewed compared to 
the adiabatic picture, to avoid crossing. In the critical range, 
the form of qa.bit) is model-independent. 



where it happens can be estimated by noticing that the 
adiabatic relaxation time t-^ (i) is a function of the in- 
stantaneous modulation phase (j) — u;pt, and (ii) sharply 
increases near the bifurcation point. As a consequence, 
sharply increases when approaches the value where qa.b 
are at their closest, because this corresponds to approach- 
ing the bifurcation point. The quasistationary (adia- 
batic) approximation requires that \dtj-/dt\ <^ 1. It is 
this condition that limits the range of adiabaticity, rather 
than a much less restrictive condition trtop <C 1. 

Even for A > A'^'^, where the adiabatic barrier in Fig.Q] 
disappears for a portion of the period, the states Qa^t may 
still coexist. A sufficiently large fluctuation is then re- 
quired to move the system away from the stable periodic 
state. It is in this region that the new scaling of the acti- 
vation energy R emerges. The control parameter is now 
J] (X Ac — A, where A^ is the "true" bifurcation value of 
the modulation amplitude. 

In the limit LOptr 3> 1 the behavior near a bifurcation 
point is in some sense simpler. In this case gait) and qb(t) 
come close to each other everywhere on the cycle, not just 
for a part of the period. The motion of the system in the 
vicinity of qa,b{t) is oscillations about the state qa{t) with 
a slowly varying amplitude. The amplitude change can 
be described by averaging the complete dynamics over 
the period. It is then mapped onto motion in an effec- 
tively stationary potential. Not surprisingly, the scaling 
of the escape activation energy R with the distance to the 
bifurcation point is the same as for stationary systems. 

In Sec. II and Appendix A we provide a general for- 
mulation of the problem of activated escape in period- 
ically modulated systems driven by Gaussian noise. In 
Sec. Ill A and Appendix B we discuss the dynamics near 
a bifurcation point in the adiabatic limit ujpty — > 0. In 
Sec. Ill B we consider the strongly nonadiabatic dynam- 



ics that emerges where still uptr <C 1. In Sec. Ill C and 
Appendix C the dynamics near a bifurcation point is de- 
scribed in the parameter range where the field becomes 
effectively fast-oscillating, i.e. uiptr ^ 1, even though the 
relaxation time in the absence of modulation ir°^ may be 
< I/llif- The connection between the nonadiabatic lo- 
cal theory of Sec. Ill B and the theory of Sec. Ill C is 
discussed in Sec. Ill D. In Sec. IV the activation energy 
is explicitly evaluated in the three regions discussed in 
Sec. Ill, and the scaling laws for the activation energy 
R (X (Ac — A)^ in these regions are obtained. The scal- 
ing crossovers are discussed. We also find nonadiabatic 
corrections to the escape rate in the adiabatic region. In 
Sec. V we consider a periodically modulated Brownian 
particle. Numerical results for the activation energy are 
compared to the results of Monte Carlo simulations and 
to the predictions of Sec. IV. Sec. VI contains concluding 
remarks. 



II. ACTIVATED ESCAPE: GENERAL 
FORMULATION 

We will adopt a phenomenological approach in which 
a multidimensional system with dynamical variables q(t) 
is described by the Langevin equation 

q = K(q;A,i)+f(<), K(q; A, i+r^-) = K(q; A, i)- (1) 

The function K is periodic in time, with the modulation 
period tf = ^tt/lof; A is a, control parameter that char- 
acterizes the modulation strength. For example, in the 
case of an overdamped particle in a potential Uq (q) mod- 
ulated by an additive periodic force F{t), the vector K 
becomes 

K(q;A,i) = -VC/o(q)+F(t). (2) 

(here and below, V = d/dq). In this case A — max|F| 
is the modulation amplitude (note that the force F{t) — 
F{t + tf) does not have to be sinusoidal). 

The function f{t) in Eq. ^ is zero- mean Gaussian 
noise with correlation matrix 

^,,it^t')^{Mt)f,{t')). (3) 

The characteristic noise intensity D can be defined as the 
maximal value of the power spectrum, 

£> = max$„„(w), $„™(w) = y die'"V«m(0- (4) 

For thermal fluctuations D — 2fcsT. The noise intensity 
D is the smallest parameter of the theory. Smallness of D 
leads to the rate of noise-induced escape W being much 
smaller than t^^ and luf- 

In the absence of noise, Eq. may have different pe- 
riodic solutions Qper, which can be stable (attractors), 
unstable (repellers), or hyperbolic (saddles). We are in- 
terested in the parameter range where one of the stable 



3 



periodic solutions q_a{t) = q_a{t + Tp) conies close to a 
saddle-type periodic solution (t) with the same period 
(period 1, for concreteness). For slow modulation, these 
states are sketched in Fig.^ They merge together at the 
saddle-node bifurcation point A = Ac. In what follows 
we will assume that A is close to the critical value Ac. 

Escape from a metastable state (t) occurs as a result 
of a large fluctuation. The fluctuational force i{t) has to 
overcome the restoring force K and drive the system away 
from the basin of attraction to qa(i) [e.g., away from the 
potential well in Fig.^. We will assume that the required 
force f (t) is much larger than the typical noise amplitude 
(X D^/\ 

The motion of the system during escape is random. 
However, different trajectories have exponentially differ- 
ent probabilities. The system is most likely to move along 
a particular trajectory called the optimal path qopt(0 
|24| . It is determined by the most probable noise real- 
ization fopt(i). In the case of a periodically modulated 
ID system driven by stationary Gaussian noise, a way 
to find the optimal paths was discussed earlier [ig. We 
now briefly outline a generalization of the formulation 
to multidimensional systems, following the arguments in 
Ref. (more details are provided in Appendix A) . 

For a stationary Gaussian noise, the probability den- 
sity of realizations of f (t) is given by the functional (cf. 
Ref. m 

(5) 



fopt(i) for t — > ±oo 



(8) 



P[m]=exp{-TZo[iit)]/D), 
where TZq is quadratic in f , 

The matrix J-' is the inverse of ipij{t — t')/D, 



dt'T^j{t - t')ipjk{t' - t") = DSikSit - t"). 



(6) 



(7) 



We are interested in noise realizations that lead to es- 
cape, and therefore f(i) largely exceeds its root- mean- 
square value. From Eq. ((SJ, the probabilities of such noise 
realizations are exponentially small and exponentially 
strongly depend on the form of f [t) . As a consequence, 
escape trajectories should form a narrow "tube" centered 
at an optimal path fopt(i) that maximizes V[f{i)\, i.e., 
minimizes 7?.o[f(t)]. The minimum of 7?,o should be found 
with the constraints that (i) the system and noise trajec- 
tories, qopt (t) and fopt {t) , are interrelated by the equation 
of motion (ii) the path qopt(0 starts in the vicinity of 
the stable state qa(i) and ends behind or on the bound- 
ary of the basin of attraction to qa [t] , and (iii) the force 
fopt(i) is equal to zero before the escape event happens 
and becomes equal to zero once the system has escaped, 
so that, as fopt(i) decays, it does not drag the system 
back to the basin of attraction to qa- 

As explained in Appendix A, these conditions lead to 
boundary conditions for optimal paths of the form 



(note that qopt(i) ends on the basin boundary, not on 
another attractor.) 

The variational problem for optimal paths is reduced 
to minimizing the functional 

7^[q,f] -7^o[f] 

+ J dt'\{t') ■ [q(i') - K(q; A,t') - i{t')] (9) 

with boundary conditions ©. The function \{t) is a La- 
grange multiplier. The boundary condition for it follows 
from the fact that the system only needs to be driven by 
a large force f it) when it moves from the attractor to the 
basin boundary. Therefore \{t) for t — *■ ±oo. 

It follows from Eqs. ©, Q and from the results of 
Appendix A that the optimal trajectories qopt (i) , fopt (i) 
are instanton-like p7| | . The typical duration of motion is 
given by the relaxation time of the system and the noise 
correlation time ^con - In stationary systems instantons 
are translationally invariant with respect to time, i.e., if 
qopt(0: fopt(i) is a solution, then qopt(< + r), fopt(i -I- r) 
is also a solution, for an arbitrary r. In contrast, in 
periodically modulated systems this is true only for r = 
Tp. The instantons are synchronized by the modulation: 
generally there is one instanton per period that would 
provide a global minimum to 7?,. 

From Eq. (O, we obtain for the escape probability 



W (xexpi-R/D), i? = min7^[q,f]. 



(10) 



qopt(i) 



qait) for t 
qb{t) for t - 



oo. 



The activation energy R is equal to the functional 
7?.o[fopt] calculated for the optimal noise trajectory for 
escape. 

For small noise intensity, the escape rate W ^ uip. It 
periodically depends on time. However, in the small-D 
limit this dependence is seen only in the prefactor [Tsl 
IT^ IT^ . Here we are interested in the exponent, which 
gives the period-averaged escape rate W. It is equal to 
the probability of escape over the time Tp divided by Tp 

In the general case, the variational problem for the ac- 
tivation energy can be solved only numerically. Therefore 
it is particularly important to find model-independent 
properties of R. So far they have been found for compar- 
atively weak modulation, where it was shown that R has 
a term linear in the modulation amplitude [l^ . In this 
paper we analyze the activation energy i? in a previously 
unexplored region near a bifurcation point and show that 
R displays a nontrivial scaling behavior in this region. 



III. DYNAMICS NEAR A BIFURCATION 
POINT 

The dynamics near a saddle-node bifurcation point has 
universal features related to the occurrence of a slow vari- 
able, or a "soft mode" For periodically modulated 
systems, closeness to the bifurcation point in the parame- 
ter space usually implies that the merging states are close 
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to each other in phase space throughout the modulation 
period. 

If the modulation frequency cup is small compared to 
the reciprocal relaxation time in the absence of modu- 
lation l/ti^\ there emerges a situation where the stable 
and unstable states come close to each other only for a 
portion of a period. During this time, the system be- 
haves as if it were close to a true bifurcation point. Then 
it is possible to single out a slow variable that controls 
the system dynamics. Escape from a metastable state is 
most probable when q^a{t) and qb{t) are closest to each 
other. 

On the other hand, if the modulation frequency ujf ^ 

l/ir°\ near the bifurcation point the states qa(i) and 
qb(t) are close to each other throughout the modulation 
period. Then escape events should be no longer synchro- 
nized by the modulation. 

Because the dynamics near a bifurcation point is slow, 
the system filters out high-frequency components of the 
noise. As a result, the noise becomes effectively S- 
correlated (we will not consider here the situation where 
the noise power spectrum has singular features at high 
frequencies). The reduction to one slow variable driven 
by white noise can be done directly in the equations of 
motion. For slow modulation {iopti^^ <^ 1), this reduc- 
tion is local in time (see Appendix B), otherwise it has 
to be done globally over the cycle (see Appendix C). Al- 
ternatively, the dimensionality reduction can be done di- 
rectly in the variational problem for the optimal escape 
path (see Appendix A). 

A. The adiabatic approximation 

In the limit of slow modulation, where the period of 
the field tf is large compared to the system relaxation 
time tr°^ , a convenient starting point of the analysis is the 
adiabatic approximation. The adiabatic periodic states 
of the system qj^^^ are given by the equation 

K(q^t;AO = o, (11) 

which is obtained by disregarding q and the noise f in 
the equation of motion ^ . 

The adiabatic stable state (attractor) qj^"^ = cia'^{t) is 
the solution qp^j. for which the real parts of the eigenval- 
ues of the matrix /t = {dKi/dqj) are all negative. These 
eigenvalues give the "instantaneous" relaxation rates, for 
a given phase of the modulation (j) = LOpt. For the peri- 
odic adiabatic saddle- type state q^'^(i) one of the eigen- 
values of pL has a positive real part. 

In the adiabatic approximation, the saddle-node bifur- 
cation occurs in the following way. At the critical value 
of the control parameter A = yl^"^, the periodic trajecto- 
ries c^{t) and c^'^{t) given by ((TTll merge, but it hap- 
pens only once per period. One can picture it by looking 
at Fig. n (b) and imagining that the states <\a{t) and 
<\b{t) touch each other. We set the corresponding instant 



of time equal to i = (or < = titf), i.e., we assume 
that q^'^(O) = q^'^(O) for Af. Additionally, we set 
qa'(O) = qf (0) = 0. 

At the adiabatic bifurcation point A — A^'^ , t = 
one of the eigenvalues /ii of the matrix pL is equal to 
zero, whereas all other eigenvalues /ii>i have negative 
real parts. The adiabatic approximation means that 
Re ^j^^i>l 2> top, or equivalently, that the relaxation time 
max|Re/Xi>i|~^ is small compared to rj?. This relaxation 

time is typically of the order of tf^ . 

We now write the dynamical variables q in the basis of 
the right eigenvectors of the matrix /t at the bifurcation 
point and expand K in the equation of motion |^ near 
the bifurcation point q = 0, i = 0, A = A^^. As shown in 
Appendix B, the motion described by the variable qi is 
much slower than the motion described by the variables 
qiyi- Over the time ~ tf'' they "adjust" to the value of 
qi, i.e., they follow qi adiabatically. The variable qi is 
the soft mode. It satisfies the equation of motion 

qi^K{q^-A,t) + h(t), (12) 
K = aql +P5A''^ - a-f'^{ujFtf. 

Here SA'^'^ = A — A^'^] the parameters a,/3,7 are ex- 
pressed in terms of the derivatives of K at the adiabatic 
bifurcation point and are given by Eqs. HB3p . (|B4p . 

The stable and unstable adiabatic periodic states in 
the absence of noise exist for a (3 5A^'^ < 0. For concrete- 
ness and without loss of generality we set a > 0. For 
small \u!Ft\ the adiabatic states can be found by setting 
K((gi):d;A,t) =0. This gives 

{q^):%=T{2atfr\ 

where t'^'^ is the instantaneous adiabatic relaxation time. 
It is given by \dK/dqi\^^ evaluated for qi = {qi)a^, 

tf ^ i [~ap6A^'^ + {ajUFtf] (13) 

[the exphcit expression for {qi)a%{t) is given by Eq. 1B5|I ]. 

The applicability of the adiabatic approximation re- 
quires not only that ajpt'^'^ < 1, but also \dt'^'^/dt\ < 
1, otherwise the system cannot follow the modulation. 
From Eq. H13|l . near the bifurcation point the time de- 
pendence of t'^'^ is pronounced, so that 

max\dtf/dt\ = 3-^/^jujf/\PSA'"^\ > ajptf. 

Therefore the inequality \dtf/dt\ ^ 1 is much stronger 
than ajptf < 1. It holds if 

tf<t:ti, ti^iajLUp)-'/^, (14) 

i.e.. Up ^ IPSA'^'^l/j. The time ti sets a new dynami- 
cal time scale that restricts the range of validity of the 
adiabatic approximation. It imposes an upper bound on 
the relaxation time of a periodically driven system where 
this approximation still applies. 
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B. Locally nonadiabatic regime 

As A approaches A^'^, the criterion H14|l is violated. 
The periodic stable and unstable states are pressed 
against each other. Since they cannot cross, they become 
distorted, as shown in Fig. Ultimately they merge, 
but along a line rather than at a point. From Eq. (|12|) one 
can see that this is just a straight line, which is described 
by the equation 

qic{t) = "fLupt. (15) 

Eqs. H12|l . IjlSI) define the new nonadiabatic bifurca- 
tional value of the modulation amplitude for slow driv- 
ing, 

Af ^ Af + p-^-fUF. (16) 

The corrections to Af of higher order in ojp are discussed 
in Sec. Ill C and obtained explicitly in Sec. V. 

It is convenient to change in Eq. H12|) to dimensionless 
variables Q — atiqi and r = t/ti, and to introduce the 
control parameter 

r, = Piju;F)-\A^^ - A). (17) 

This is the only parameter of a slowly driven system. It 
describes both adiabatic and nonadiabatic behavior and 
gives the reduced distance to the bifurcation point. 

The equation of motion for the reduced variable takes 
the form 

^ = G(Q, 77, r) + /(r), G = - + 1 - ,7. (18) 

The function /(r) = (7Wf)"Vi(0 in Eq. (O de- 
scribes reduced noise. It is effectively 5-correlated on 
the slow-time scale, (/(t)/(t')) — 2D6{t — r'), as ex- 
plained in Appendix B. The noise intensity D is given by 
Eq. iBTll . 

The stable and unstable states Qa,b{T) are given by 
the equation dQ/dr = G. In the absence of noise this 
equation has symmetry Q ~Q,t — r. As a conse- 
quence, the stable and unstable states are antisymmet- 
ric, (5b(r) = —Qa{—T). Therefore it suffices to find only 

Qa{T). 

We start with the adiabatic approximation. It applies 
for 7] ^ 1. The adiabatic stable and unstable states in 
the reduced variables are given by the equation G = 
and have the form 

Qtb = T[r' + {^-lf]"'. (19) 

Each of these states is symmetric with respect to r = 
0, where they are closest to each other. The adiabatic 
bifurcation point is 77 = 1, which corresponds to A = A^*^. 

The region 77 1 is nonadiabatic, and 77 = (or 
A = j4^') is the nonadiabatic bifurcation point for slow 
driving. At this point Qa,b{T) merge into the straight 
line Qc{t) = t. 



Close to the nonadiabatic bifurcation point, where 
7; <C 1, one can find Qa,b{'T) by perturbation theory in 
the whole range —00 < t < |ln7;|^/^/2. The linearized 
equation for the difference 6Q(t) = Qa{T) — t has a sim- 
ple form = 2t SQ + 77. By solving it we obtain 

Qa{T) = -Qb{~T) « r - 77/:^ dn e-'-?. (20) 

In the region of large negative r the function Qa (r) = 
—Qb{~T) has a simple form Qa{T) « — r — 77(2t)~^. The 
states Qa and Qb are closest to each other, with separa- 
tion ^ 77, in the range |t| < |ln7/|^/^/2. The interstate 
separation decreases as rj approaches the bifurcational 
value rj = 0. At the same time, the range of r where Qa 
and Qb are close to each other increases with decreasing 
V- 

As T increases beyond w | In 771^/^/2, there occurs a 
sharp crossover from the nearly linear in r solution for 
Qq(t) (I^OJ) to the adiabatic solution Qa oc -r. The 
functions Qa,b{T) for a specific value of 77 are shown in 
Fig. El 
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FIG. 2: Nonadiabatic stable and unstable states Qa(r) and 
Qb(j) ~ -~Qa{—T) for slow modulation as given by the equa- 
tion dQ/dr = G{Q,r],T) for rj = 0.2. The functions Qa,b(T) 
are strongly asymmetric, in contrast to the adiabatic states 
I19II which are even functions of r. 

The interval of the real time \t\ ^ tr = t;|ln7;|^/^, 
where the states Qa,b are nearly linear in t, should be 
much smaller than I/ojf in order for Eq. H12|l to apply. 
This imposes a restriction on rj, 

77 > exp(-G|a|7/wj.), G-1. (21) 

For smaller 77 cx \A^^ — A\ the local approximation, where 
the coefficients are expanded about the adiabatic bifur- 
cation point, no longer applies. The relaxation time 
becomes comparable to the modulation period. It follows 
from Eq. (|21|) . however, that for low frequencies the local 
approximation is extremely good. 

On the whole, the locally nonadiabatic regime is lim- 
ited in rj by the condition 7y ^ 1 and by Eq. (|21|) . The 
width of the amplitude range A^' — A imposed by the first 
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condition linearly increases with the field frequency, in 
the approximation Therefore locally-nonadiabatic 

critical behavior is more pronounced for higher frequen- 
cies. However, the appropriate frequency range is limited 
from above by the condition H21() . For higher ujp there 
should occur a crossover to a fully nonlocal picture, which 
is discussed in the next section. 



C. Fast-oscillating field 

Sufficiently close to the "true" critical value of the 
modulation amplitude Ac, the relaxation time of the sys- 
tem becomes large compared to the modulation period 
even if Loptr^^ ^ 1 far from the bifurcation point. The 
inequality upt,. 3> 1 defines the third region, in addition 
to the adiabatic and locally nonadiabatic, where we have 
analyzed the dynamics near a bifurcation point. The 
analysis of this region is simplified by the fact that here 
the modulating field is effectively fast-oscillating. 

For LUptr 2> 1, near the bifurcation point the periodic 
stable and unstable states q_a{t) and qb(i) stay close to 
each other throughout the cycle, see Fig. 13 For A = Ac, 
they coalesce into a periodic critical cycle qc(i) = qc(i + 
Tp). When A is close to Ac and q is close to qc, we can 
simplify the equations of motion |^ by expanding the 
function K in (5q = q — q^ and SA — A — Ac (cf. Ref . ITlh . 

^q = /i5q+i(5q. V)2K + &4^K + f(t). (22) 

Here /i^ = (Xij{t) = dKi/dqj. All derivatives of K are 
now evaluated for A ^ Ac and q — (lc{t). Therefore all 
coefficients in Eq. ()22|l are periodic functions of time. 




FIG. 3: The stable and unstable states, qa(t) and q6(t), close 
to the bifurcation point. For Ldpti ^ 1 the states are close 
to each other throughout the modulation period. The figure 
refers to a one- dimensional overdamped particle in a potential 
U{q,t) = iq-ig^-Aqcos(wFt) for {Ac -A) /Ac ~ 0.01. The 
modulation is comparatively slow, ajF^r"' = 1, but for chosen 
A the relaxation time becomes long, uiptr ~ 9.8. 

If initially the system is close to qdt), its distance from 



qc(t) will oscillate with frequency ujf and with an ampli- 
tude that slowly varies over the period Tp. This ampli- 
tude is a slow variable, Q^'^{t). The equation for Q^'^{t) 
can be obtained by an appropriate averaging method ex- 
plained in Appendix C. After rescaling to dimensionless 
coordinate Q cx Q^'" and time t cx t, see Eq. (|C11|) . this 
equation takes a form which is similar to Eq. I|18|l . 

^-G(Q,r7) + /(r), (23) 

G^Q^-V, rj^(3'{Ac-A) 

[in contrast to Eq. H18() . the function G here is indepen- 
dent of time] . The coefficient /3' is given by Eq. HC7() . 

The parameter rj in Eq. (|23|l is the scaled distance 
to the bifurcation point. The stationary states Qa.b = 
=F77^/^ exist for 77 > 0. They merge for 77 = 0. The noise 
/(t) is effectively white on the time scale that largely 
exceeds ujp^ and the noise correlation time tcorr- Its in- 
tensity D is given by Eq. (|C12|I . 

The results of this Section and Appendix C refer to the 
case (jjpt^- 3> 1, but arbitrary ujpti^ . Therefore the prob- 
lem is different from the standard problem of slow mo- 
tion in a fast-oscillating field , where of interest is the 
smooth term in the oscillating coordinate. In contrast, 
here we are interested in the slowly varying oscillation 
amplitude. If ajpti^^ » 1, a transition to slow and fast 
variables can be made already in the original equation of 
motion ^ , by separating q into slow and fast oscillating 
parts. The equation for the slow part near the bifurcation 
point will again have the form H23() . but the expressions 
(|C7p for a',P' will be simplified [in particular, the factor 
Kii in Eq. (|C7|I will be equal to one]. 

D. Connection to the locally nonadiabatic regime 

Eq. H22|l allows us to look from a different perspec- 
tive at the locally nonadiabatic regime that emerges for 
LUpti- ^ 1. In contrast to the approach of Sec. Ill B, 
where the starting point was the adiabatic approxima- 
tion, here we will assume that A is close to the true bi- 
furcational value of the amplitude Ac and that q(f) is 
close to the critical cycle qc(t), at least for a part of the 
period Tp. 

For LUptj: <g; 1, one can think of a local in time de- 
scription of the dynamics near the cycle qc(0- From 
Eq. H22|) , this dynamics is determined by the eigenvalues 
/i^(i) of the matrix /t(t). In contrast to the analysis of 
Sec. HI A, we consider here the matrix jj calculated for 
the critical cycle qc(t) rather than the two similar ma- 
trices calculated separately for the adiabatic stable and 
unstable states. 

For much of the driving period the real parts of iiu{t) 

are all large, [Re /Lt^l ~ l/tr°^ » ojp. Then, when the 
system is in the stable state, it follows the field adia- 
batically. The adiabaticity is broken where one of the 
eigenvalues, say ^i(t), goes through zero. As we will see. 
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at this time the stable and unstable states are most close 
to each other. We set the time when it happens equal to 
zero, i.e., ^i(O) = 0. 

For small \t\ <C Tp the analysis of the system dy- 
namics is in many respects similar to that in Sec. Ill B 
and Appendix C. First, Sq{t) in Eq. H22|l is written as 
(5gi/ej^(0), where e^{Q) are the right eigenvectors of 
the matrix //(O). The component Sqi of (5q along the 
eigenvector ei(0) of /i(0) will be the slow variable, or the 
soft mode. 

The matrix fl(t) can be expanded about t = for small 



m « A(o) + /i(o)t, 



(24) 



where the time derivative is taken for t — 0. This deriva- 
tive is small, its matrix elements on the eigenvectors 
e^(0),e^,(0), are - up/i^'^ < [here, e^(0) 

are the left eigenvectors of the matrix /i(0)]. 

With (|23Il, Eq. (|22Jl can be solved for the "fast" compo- 
nents 6q^yi. Over a short time ~ ir"' they approach their 
quasistationary values for given Sqi. Those are small, of 
order Sqf,SA, Sqiujpt, and follow 6qi adiabatically. Noise- 
induced fluctuations of Sqi,^i about the quasistationary 
values are also small for small noise intensity. Therefore 
the effect of Sq^yi on the dynamics of Sqi can be disre- 
garded. 

The equation of motion for 6qi is of the form 



6qi = fj-itSqi + aSqi 
hit) = ei(0).f(t), 



h/5M + /i(i), (25) 
/ii =ei(0) -^(0)61(0). 



Here, a = (l/2)(ei(0)-V) Xi, [3^dKi/dA, with Ki = 
§1 (0) • K being now the component of K in the direction 
ei(0). All derivatives of K are calculated on the critical 

cycle qc(t) for t = 0. Because ~ ujF/t'f^ is small, 
relaxation of Sqi is slow compared to relaxation of 6qi,yi, 
for typical \t\ <^tf- 

Eq. H25|) describes the stable and unstable states of the 
original equation of motion in the region |i| ^ r^^. It 
is seen that these states, {Sqi)a and (Sqi)i,, exist provided 



/ii > 0, a/3 (£4 < 0. 



(26) 



In this range Eq. (^51) is equivalent to Eq. ((TH|l . This 
can be seen if, on the one hand, Eq. (|18|l is written for 
the deviation SQ = Q — t of Q from its value on the 
critical cycle Qc — r, and on the other hand, in Eq. H25|) 
one changes to scaled variables SQ = a{2/ jiiY^^Sqi and 
T = (/ii/2)^/^t. The control parameter r] in (|18() becomes 



-2a(3^i^5A. 



(27) 



The analysis of Eq. H18|l then applies to Eq. (|25|) . In par- 
ticular, the statement in the beginning of this subsection 
that the stable and unstable states are at their closest 
for f = is an immediate consequence of the explicit 
expression for these states (|2(JII . 



There is an important difference between this approach 
and the approach of Sec. Ill B. Because here we do not 
start from the adiabatic approximation, we formally have 
not specified how small is the difference between the criti- 
cal amplitude Ac and its adiabatic value A^'^. In Eq. H16() 
we only obtained the linear in ujpti^'^ term in Ac^A^'^. In 
general, Ac — A^'^ has also higher-order terms. They can 
be obtained by taking into account the dependence of the 
coefficients a, /?, 7 in Eq. ljB3|) on A, which was previously 
disregarded. This is illustrated for a particular model in 
Sec. V. It is for the renormalized critical amplitude, i.e., 
for the control parameter given by Eq. H27|) . that there 
holds the exponential limit (|21|l on the range where the 
local nonadiabatic approximation applies. The inequal- 
ity (|21|l indicates that, for small frequency, the critical 
amplitude found from the local theory is exponentially 
close to the exact Ac- This is confirmed by numerical 
calculations for a model discussed in Sec. V. 



IV. ACTIVATION ENERGY OF ESCAPE 

It follows from the results of Sec. Ill that, near a bifur- 
cation point, a periodically driven system has a soft mode 
Q, and the noise that drives this mode is effectively white. 
The equation of motion is of the form dQ/dr — G + /{t) 
l(TH|l . where the function G is given by G = Q^-l-l — 77 — 
for Luptr < 1 [cf. Eq. (jTHJ)] and G = - 77 for ajpti > 1 
[cf. Eq. 123)] ■ The intensity D of the noise /(t) has the 
form (jB7p and (|C12p in these two cases. 

For a white-noise driven system, the variational prob- 
lem IjlUI) of calculating the period-averaged rate of 
activated escape W can be written in the form 

W = const X exp{-R/D), R = min7^[Q], (28) 

(cf. App endix A). In contrast to the standard formula- 
tion [29j , the function G here may depend on time r and 
is not time-periodic, in the actual range of r. The min- 
imization is carried out over the paths Q{t) that start 
at Qa{T) for r — > — 00 and end at Qb{T) for r — » +cx3, 
with time-dependent stable and unstable states Qa,b{T)- 
The non-stationarity emerges for slow modulation, where 
LUFt-c ^ 1, and is related to the assumptions that (i) es- 
cape is most likely to occur during a portion of the period 
where the states Qa.b are close, and (ii) the duration of 
motion along the optimal escape path Qopt{T) is much 
less than the modulation period. 

We have solved the variational problem using the 
Hamilton equations of motion for Q and P = 
dL/d{dQ/dT), 

We then verified the assumptions made in obtaining 
Eqs. (EHIl, 
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Equations (|29|l were solved numerically. This was done 
by choosing the initial conditions on the optimal path 
close to Qq(t) with large but finite negative r. In this 
range Eqs. (|29|l can be linearized in Q — Qa{T). On the 
solution that goes away from Qa the momentum P is 
linear in Q — Qa, 

P^[Q-Qa{r)]/al{T) (30) 

a2(r) = 2 / dr' exp 4 / dr" Qair") 

J — oo \_ J t' 

We used the shooting method: we sought such initial 
Q — Qa{T) for given r that the trajectory approaches 
Ob(T) for large r, cf. RefEl 

Numerical results for the activation energy in the whole 
range of slow driving, where LOptv 1, are shown below 
in Figs.0IIHon linear and logarithmic scales, respectively. 
Note that the activation energy is a function of a single 
control parameter 77 oc Ac — A, and in this sense the 
results are universal, i.e. system-independent. In the rest 
of this Section we discuss analytical results and compare 
them with the numerical results. 



A. Activation energy in the adiabatic 
approximation 

The adiabatic regime applies when the driving is slow, 
ti^F^r ^ 1, and the system is sufficiently far from the 
bifurcation point, so that ^ ti [cf. Eq. (|14|l ] or equiv- 
alently |(5A^^| = \A'f - A| > lof\iIP\- Respectively, the 
dimensionless control parameter rj (x Ac — A is large, 
?7 — 1 3> 1, [cf. Eq. H17(l : we note that the actual parame- 
ter in the adiabatic range is not rj but ?7 — 1] . In this case 
we expect that escape occurs when the adiabatic states 
(|19|l are closest to each other, i.e., for r = 0. Then, in 
the first approximation, the term in the function G in 
Eq. H18|) can be disregarded, and G becomes 



(31) 



The solution of Eq. with G of the form and 
with boundary conditions Q{t) — > ^{rj — 1)^/^ for r 
=poo is well known. It is an instanton (kink) 



QlU^, ro) = {V- 1)'/' tanh[(7y - l^/^r To)] 



(32) 



centered at an arbitrary tq. 

The characteristic duration of motion along the path 
QoptC''') '^^ dimensionless time is At ^ {rj — 1)^^/-^, which 
corresponds to At ~ t-c in real time. Since At <^ {rj — 
1)^/^, the term t^ in the function G [Eq. (^HJl] can be 
disregarded compared to 77 — 1, which justifies replacing 
G with G^'i as long as |to| < (?? - 1)^/^ 

The activation energy (|28|l calculated along the path 
Qopt is 



This equation shows that the activation energy of escape 
scales with the distance to the bifurcation point as {Ac — 
A)^ with ^ = 3/2 in the adiabatic region. 



B. Nonadiabatic correction to the activation 
energy 

We now consider the lowest-order correction to the adi- 
abatic activation energy. Two factors have to be taken 
into account. First is that, because of the nonzero dura- 
tion of motion along the escape path At, the equilibrium 
states Q'^'^i,{t) change, which was disregarded in the anal- 
ysis of Sec. IV A. However, the corresponding correction 
to R is exponentially small. Indeed, if we consider TZ as 
a function of the end point Q on the optimal path, we 
have = |P|, where P is the momentum on the 

optimal path. For an instantonic solution, the momen- 
tum goes to zero as Qopt Qa.b, see Eq. (fSnil . Therefore 
a small change of Qa.t affects the activation energy very 
weakly. 

The major nonadiabatic correction to R comes from 
the time-dependent term in G = G'''^ - t^ [cf. Eq. ((TH|) ]. 
This term lifts the time invariance of the instanton 
QoptiT,TQ) with respect to tq. 

To first order in t^, i.e., to lowest order in (ry— 1)~^, the 
correction 5R can be found from Eq. H28|l by integrating 
the term t^ along the zeroth order trajectory Qopt(''': '^0), 



6R 



mm / dr t . 

To J dr' 



(34) 



(we used that dQ^'^^/dr = — G^'*). Minimization here is 
done over tq , the center of the instanton. It is necessary 
because R is the absolute minimum of the functional TZ. 

A direct calculation shows that the minimum of 5R is 
reached for tq = 0, and 



5R^^-{^-l)-y\ 



(35) 



R 



■ad 



(r,-l)3/2cx(Af -A) 



\3/2 



(33) 



The correction 5R rapidly falls off with increasing rj — 1. 
On the other hand, as 77 decreases and becomes ^ 1, the 
term SR increases very fast, which indicates a breakdown 
of the adiabatic theory in this region. 

The analytical results in the adiabatic region are com- 
pared with the numerical solution for the activation en- 
ergy R in Fig. ^ The corrected adiabatic theory works 
well in the whole range where the control parameter 
rj ^ 3, but for smaller 77 nonadiabatic effects are sig- 
nificant and have to be taken into account in a nonper- 
turbative way. 



C. Activation energy in the locally nonadiabatic 
region 

Standard techniques do not allow to solve equations 
(|29|l for the optimal path analytically in the general case 



FIG. 4: The activation energy R vs. rj cc Ac — A for slow 
driving, ujFtv <S 1- The thick solid line shows the numerical 
solution of Eq. I|28|l . The dashed line is the adiabatic activa- 
tion energy (13^ . R'"^ (x {v — 1)^''^. The thin solid line shows 
the corrected adiabatic activation energy R'^'^ + SR. It is close 
to the numerical result for rj ^3. The correction SR diverges 
at the adiabatic bifurcation point rj — 1. 

r; ~ 1. This is because the function G in Eq. 129|l ex- 
plicitly depends on time r. However, a solution can be 
obtained close to the bifurcation point, where rj oc Ac — A 
is small, but not exponentially small, cf. Eq. (|21|l . 

Unusually for an instanton-type problem, and because 
of the strong time dependence of G, the optimal path 
can be found by linearizing the equations of motion H29() 
about the critical state Qc = t. This gives 

SQ ^ 2t6Q -r] + 2P, P = -2Pt, (36) 

where 6Q = Q — t. The solution of these equations with 
boundary conditions Q{t) Qa,b{T) for t =Foo is 

Oopt(T) = T-T] f dT'[l-V2e-^'"]e^"-^'\ 
Jo 

Popt{T)=ve-^'/V2, (37) 
where we took into account the explicit form of Qa.b{t) 

It is seen from Eq. 1)3 7|) that the momentum on the 
optimal path Popt has a shape of a Gaussian pulse cen- 
tered at r = 0, with width ~ 1. The coordinate Qopt(''") 
over the dimensionless time r ^ 1 switches between the 
equilibrium states Qa,b- The typical duration of motion 
in real time is ti . 

From Eqs. (|37(l . the nonadiabatic activation energy of 
escape for w^ir <C 1 is 

^nonad ^ (7^/8)1/2,^2 ^ _ ^^2 _ (gg) 

Here, the critical amplitude Ac is given by Eq. to 
first order in ujp. 

It is seen from Eq. H38() that, in the locally nonadia- 
batic region, the activation energy again displays scaling 
behavior, R oc {Ac — A)^. But the scaling exponent is 
5 = 2, it differs from the adiabatic exponent ^ = 3/2 



FIG. 5: The activation energy R — — D In VK on a logarithmic 
and linear scale (inset) vs. ij cc Ac — A for slow modulation, 
Lvpti <^ 1. Thick solid lines show the numerical solution of 
the variational problem 12811 . It describes the crossover be- 
tween different scaling regions. The thin solid line shows the 
adiabatic scaling for large rj, R cc rf^"^. The full result of the 
adiabatic approximation is shown by the dashed line. The 
dash-dot line shows the nonadiabatic result 13811 that applies 
for »7 ^ 1; here R cc rf' . 

(|33|l that has been known for stationary systems. This 
is a result of the complicated nonadiabatic dynamics as- 
sociated with avoided crossing of the stable and unstable 
states. The onset of this scaling behavior is the central 
result of the paper. 

The predicted ^ = 2 scaling is compared with the result 
of the numerical calculation in Fig. |31 The analytical 
and numerical results are in quantitative agreement in 
the whole range rj ^2. 



D. Activation energy for Loptr 2> 1 

It was shown in Sec. HI C that, sufficiently close to a 
bifurcation point, there holds a condition uopt-c 3> 1, even 
where the modulation frequency is less than the relax- 
ation rate far away from the bifurcation point, LDpt^^^ ^ 
1. Finding the activation energy of escape R = R^^'^ for 
LOpti: 3> 1 is formally similar to that in the adiabatic ap- 
proximation. The only difference is that rj—1 in Eq. H31fl 
should be now replaced by rj. This gives 

^fa.t ^ (4/3)^3/2 ^ _ ^^3/2^ (39) 

Both the coefficient /3' that relates rj to Ac — A [see 
Eqs. (|23l and HC7|l ] and the noise intensity D HC12|I 
depend on the arbitrary initial time tj. It enters the 
weighting factor Kii(i,ti) which was used in obtaining 
the equation of motion for the slow variable H23|l . A 
straightforward analysis shows that ti drops out from 
the ratio l3'^/'^/D, which gives the escape rate W oc 
exp(-^f''''Vi)). 

Eq. H39|) shows that the activation energy displays scal- 
ing behavior with the distance to the bifurcation point 
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in the range LOptr ^ 1- The scahng exponent is f = 3/2, 
as in the adiabatic case. 



E. Scaling crossovers near a critical point 

Eqs. (|33|l . H38|l . and (|39|l show the onset of three re- 
gions where the activation energy of escape displays scal- 
ing dependence on the modulation amplitude, R oc R oc 
{Ac — A)^. The adiabatic and locally nonadiabatic re- 
gions emerge only if the modulation frequency is slow 
compared to the relaxation rate far from the bifurcation 
point, ajptf'^ 1- In this case, as seen from Fig. |S1 as 
the bifurcation point is approached, the system displays 
first the adiabatic scaling ^ = 3/2, which for smaller 
Ac — A goes over into the scaling ^ = 2. As the bifur- 
cation point Ac is approached even closer, there emerges 
the fast-oscillating regime where ^ = 3/2 again. 

The widths of the regions of different scaling strongly 
depend on the modulation frequency. For small tupti^^ «C 
1 the range of amplitudes where motion is effectively fast 
oscillating, wj^ir 3> 1, is exponentially narrow. However 
this range increases very rapidly with increasing ujp. The 
particular way in which the widths of different scaling 
regions vary with uf depends on the system dynamics. 
Ultimately, for ujptf'^ > 1, the regime of effectively fast 
oscillations becomes the only observable regime near a 
bifurcation point. 



bifurcational value of the coordinate is Qc = 0, and the 
adiabatic bifurcation occurs for t = (or equivalently, 
t — riTp). Near the adiabatic bifurcation point we have 

q = q^ + 5A^''-]^Au^ie+f{t), (41) 

with SA^'^ — A — A^^. This equation will have the same 
form as Eq. H12|l if we replace the factor A in Aupt'^ with 
Ac'^ = 1/4 [as it was done in Eq. lO]. 

From Eq. (|41|l it follows that, for the model under con- 
sideration, the parameters in Eq. (|12|l are a — (3 = 1,7 = 
{Af/2y/^ = 2-3/2. Therefore, from Eq. lUnjl, to first or- 
der in ujp the critical amplitude is Af = 1/4-1- 2~^^^u;f. 
It is shown in the main part of Fig. |^ with the dashed 
line. 




V. SCALING CROSSOVERS FOR A MODEL 
SYSTEM 

To test the occurrence of three scaling regions and the 
scaling crossovers, we have studied activated escape for 
a model system, an overdamped Brownian particle in a 
modulated potential well. It is described by the Langevin 
equation 

. ^ „5C^ ^ ^^^^^ if it) fit')) = 2D5{t - t'), 



dq 

U{q,t) = ^g?"^ ^ 4'^ ~ AqcosLUpt. 



(40) 



The shape of the potential U {q, t) is shown schemat- 
ically in Fig. n In the absence of modulation, A — 0, 
the system has a metastable state at the bottom of the 
potential well, qa = ^1/2, and an unstable equilibrium 
at the barrier top, qb = 1/2. The relaxation time is 
4°^ = 1/U"{qa) = 1. In the presence of modulation, the 
states qa.b oscillate in time. As the modulation ampli- 
tude A increases to the critical value Ac (the saddle-node 
bifurcation), the states merge, and then, for A > Ac, dis- 
appear. 

The frequency dependence of the critical amplitude Ac 
is shown in Fig. In the limit ojp = we have Ac = 
Af = 1/4. The hnear in LOp correction to Ac can be 
obtained from Eq. (|16|) by noticing that the adiabatic 



FIG. 6: The critical amplitude Ac as a function of the mod- 
ulation frequency luf for the system (1401 . Numerical results 
are shown by thick solid lines. The dashed line shows the lin- 
ear in ojf nonadiabatic correction to Ac described by Eq. 1161 . 
The thin solid line in the inset describes a correction obtained 
from the self-consistent local analysis, Eq. 1)42^ . The dash-dot 
line describes the high-frequency asymptotic that follows from 
Eq. 

As discussed in Sec. Ill D, the local nonadiabatic the- 
ory allows us to find higher-order terms in the critical 
amplitude. This is done by noticing that the critical state 
qdt) into which the stable and unstable states merge 
at the bifurcation point is linear in t for small i, i.e., 
qcit) — qc'^ cx t. By substituting this solution into Eq. H41() 
(without noise) we obtain 



Ac 



1 + u}f+ujf{cjI + 2Y^^ /A,ujpU^l. (42) 



This equation is in good agreement with the numerical 
data for ajpt^'^'' ^ 0.4, as seen from the inset in Fig. 
The difference between the numerical and analytical Ac 
decreases exponentially fast with decreasing up. 

We also evaluated for slow driving the time derivative 
of the eigenvalue /^i = 2qc{t) on the critical cycle qdt). 
For the model H41|l the stable and unstable states are at 
their closest for t = 0. Eqs. ^7^, ((23) show that, at 
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this time, /ii = {2Ac)^/^u}f- This value agrees with the 
numerical values of /ii to better than 10% for cup < 0.5. 

In the high frequency limit, WF^r"'' 3> 1, the motion 
of the system (|40|l is a superposition of slow motion and 
fast oscillations at frequency tup- To lowest order in ujp^ 
we have q ~ Q + (A/ujp) sinujpt. The equation for the 
slow variable Q becomes 

Q = Q'-l + ^ + m- (43) 

It shows that, for large ujp, we have Ac ~ Lijp/^/2. This 
is in good agreement with numerical data in Fig. IHl for 

In the intermediate range, ujpt'f^ < 1, the motion may 
not be separated into slow and fast-oscillating for weak 
modulation, but the separation becomes possible near a 
critical point, uipt^ 3> 1. Here, the coefficients in the 
equation of motion for the slow variable and the effective 
noise intensity ljC7|l . I|C12|) are nonlocal and had to be 
evaluated numerically as functions of wf- 

A. Activation energy 

For a periodically modulated overdamped Brownian 
particle described by Eq. (|40|l . the activation energy of 
escape R can be found from the variational problem lO, 
(fTn|l or, equivalently, (|2H1) ^18]. The variables Q,P,t, 
and the function G in the Lagrangian L lEHl and the 
Hamilton equations should be changed to q,p, t, and 
~dU{q, t)/dq, respectively. As explained in Sec. II, there 
is one optimal path per modulation period. The initial 
condition for the momentum p on the optimal path is 



It is seen from Fig. [7| that, for small wf, the adia- 
batic approximation applies over a broad region of driv- 
ing amplitudes. Near the bifurcation point it gives scaling 
R cx {Ac-Af with C = 3/2 (cf. the panel for ujp = 0.1). 
However, close to the bifurcation point this scaling does 
not work. Instead there emerges the nonadiabatic dy- 
namic scaling with ^ = 2. For ujp = 0.1 the range of the 
nonadiabatic scaling is comparatively narrow. 

As the frequency increases, the amplitude range char- 
acterized by the ^ = 2-scaling dramatically increases. For 
ujp = 0.25 this is practically the only scaling seen near 
the bifurcation point. 

With further increase oi up, close to the bifurcation 
point there emerges a region of the fast-oscillation scaling 
R cx {Ac — A)^ where again f — 3/2. The panel for 
ujp = 0.5 shows a crossover from the scaling ^ = 3/2 very 
close to the bifurcation point to the scaling ^ = 2 further 
away from Ac. Note that the frequency cop — 0.5 is 
neither small nor large, and therefore there is a noticeable 
difference in the coefficients at {Ac — A)'^ obtained from 



given by Eq. (|^ . with Qa{T) replaced by qa{t) [the ex- 
pression for (T^ can be further simplified taking into ac- 
count the periodicity of qa{t)]- Then Eqs. (|29l) can be 
solved numerically. 

The obtained activation energy R as function of the 
modulation amplitude for four characteristic values of 
wf is shown in Fig. The solid lines on this plot cor- 
respond to the results of the numerical solution of the 
variational problem. The dashed lines in the panels for 
Up = 0.1,0.25 (we remind that 4°'' = 1) show the adia- 
batic approximation, 

R^^ = min [U{q,{t)) C/(g,(t))] = ^ - . 

The dash-dot lines on all panels show the locally nona- 
diabatic approximation near the bifurcation point, which 
gives 

1/2 / \ "'"/'^ 

In plotting this expression we used the values of Ac and 
/ii found numerically (they were very close to the ana- 
lytical expressions given above). 

Finally, the dashed lines in Fig. [7| in the panels for 
Up = 0.5, 1 show scaling for the effectively fast-oscillating 
regime near the bifurcation point, with 

3 D 

The coefficients /3' and D as given by Eqs. (|n7|) . l(nT2|) 
were obtained numerically. 

I 

the full variational problem for R and from the locally 
nonadiabatic theory near the bifurcation point. 

When upt? ^ 1 we do not expect to see scaling for 
either adiabatic or locally nonadiabatic regime. The only 
scaling to be expected near the bifurcation point is the 
fast-oscillation one, with ^ = 3/2. This is indeed seen 
in the panel for wf — 1 in Fig- 13 We note that the 
global fast-oscillation approximation H43|) does not apply 

for upt? = 1. 



B. Simulations 

An additional test of the results can be obtained by 
directly simulating the Brownian dynamics described by 
Eq. H40|l . We conducted such simulations using the sec- 
ond order integration scheme of stochastic differential 
equations developed in Ref. H^- As a result of the sim- 
ulations we obtained the probability distribution of the 
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FIG. 7: The activation energy of escape R vs. modulation amplitude A on the logarithmic and linear (inset) scales for a 
Brownian particle in a modulated potential 14U1 . The values of ljf are indicated on each panel. The thick solid lines show 
the results of the numerical solution of the variational problem for R. The dashed lines for lof = 0.1,0.25 show the adiabatic 
approximation, whereas for ljf ~ 0.5, 1.0 they show the approximation of effectively fast oscillations: in the both cases the 
scaling exponent is ^ = 3/2 (for ljf = 0.1 this asymptotic scaling is shown by the thin solid line). The dash-dot hnes show the 
^ = 2 scaling 13811 . The dots show the results of numerical simulations of Eq. 14011 . 



dwell time in the metastable state pdw{t)- It gives the 
probabiUty density (over time) for a system prepared at 
t = close to the attractor to stay in the basin of attrac- 
tion until time t and leave at that time. 

In practice we calculated Pdw{t) by detecting the sys- 
tem at time t at a point q that lied well beyond, but not 
too far from, the oscillating boundary qb{t). It is simply 
related to the time-dependent escape probability W{t), 
which gives the probability current from the attraction 
basin at time t if the system was in the stable state at 

i = fiEl, 



Pdw(i) = W^(i)exp 



- / dt'Wit') 



(44) 



The average escape rate is given by the mean dwell time, 

-1 



W 



dttpdw{t) 







(45) 



We studied small noise intensities so that W <^ ojp', then 
W was independent of the position of the "observation" 
point q. 

In most simulations the system was prepared initially 
at the stable state qa{t); we found that W was indepen- 
dent of the initial state provided it was close to qa{t)- 

The dwell-time distribution for a particular set of pa- 
rameter values in Eq. H40I) is shown in Fig. |S| The data 
refer to modulation at a comparatively low frequency and 
with comparatively small amplitude. The function Pdw(t) 
is strongly modulated in time, with period Tp. This 
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FIG. 8: The scaled probability density of the dwell time 
Pdw{t) = Pd^{t)TF obtained by numerical simulations of a 
Brownian particle in a modulated potential, Eq. 1401 . The 
parameters are A = 0.1, D = 0.05, ojf = 0.25 (solid line). 
The dashed line shows the exponential fit of the envelope 
with decrement Wtf — 0.008. 



means that escape events are strongly synchronized by 
the modulation, in agreement with the analytical results 
for W{t) obtained for the same model in Ref. 

The average escape rate W was found from the data of 
the type shown in Fig. [S] by calculating the mean dwell 
time (|45f) and also from an exponential fit of the envelope 
of Pdw(^)- These two approaches gave the same result. 
For each set of A, ujf, D we observed ^ 10^ escape events. 
Then D was changed. The activation energy of escape 
was found from In W for 2-4 values of D. We tested that 
it was independent of D in the range R/D > 6. 

The data of simulations are shown in Fig. Qby dots. 
For all parameter values they are in excellent agreement 
with the results of the numerical solution of the varia- 
tional problem 0. 



VI. CONCLUSIONS 

In conclusion, we have identified three regions near a 
bifurcation point where the activation energy of escape 
displays scaling behavior as a function of the amplitude 
of periodic modulation. The main results refer to slow 
modulation, where LUptj- <^ 1. We show the emergence of 
nonadiabatic behavior in this region. The nonadiabatic- 
ity leads to a crossover from the scaling with exponent 
^ = 3/2, previously found for stationary systems, to a 
new dynamical scaling with ^ = 2. The — 2 region 
first emerges near the bifurcation point and then expands 
with increasing modulation frequency. With further in- 
crease oi ujp the crossover ^ = 2 to ^ = 3/2 can be ob- 
served. Again, the effectively fast-oscillating region with 
^ = 3/2-scaling emerges first near the bifurcation point. 
Even though the widths of the regions of different scaling 
depend on the parameters of a system, the phenomenon 
of scaling crossovers should be universal. 



The onset of the ^ = 2-scaling is a consequence of 
the slowing down of motion near a bifurcation point. 
The adiabatic relaxation time of the system strongly 
depends on the distance to the bifurcation point, tj- oc 
{Ac — A)~^l'^ . The nonadiabatic scaling emerges where 

tr becomes ^ [rj;-tr°'']^/^. For smaller |Ac — ^| the depen- 
dence of ovl Ac — A becomes weak, while tr still largely 
exceeds the modulation period Tp. This is associated 
with avoided crossing of the stable and unstable states, 
which occurs with decreasing A^ — A z& these states are 
pressed against each other when the system approaches 
the adiabatic bifurcation point, see Fig.|21 

Both in the adiabatic regime and the nonadiabatic 
regime for ixipt^ <C 1 escape is most likely to occur while 
the equilibrium states are close to each other. The es- 
cape rate is therefore determined by the behavior of the 
system for a small portion of the modulation period, i.e., 
locally in time. 

The regime of effectively fast oscillations near the bi- 
furcation point emerges for LOpt^ 3> 1. It can arise even 
where the modulation period Tp exceeds the relaxation 
time far from the bifurcation point. In this regime the 
motion is controlled by a slow variable, but the dynamics 
of this variable is no longer determined by local (in time) 
characteristics. The relevant characteristics are obtained 
by averaging the appropriate parameters along the criti- 
cal periodic trajectory of the system into which the sta- 
ble and unstable periodic cycles merge at the bifurcation 
point. 

We have developed a general formulation of escape of 
periodically modulated systems driven by colored Gaus- 
sian noise. Near a bifurcation point slow motion of the 
system effectively filters out high-frequency components 
of the noise spectrum and makes the noise effectively 
white. From the theoretical point of view it is interesting 
that, in the locally-nonadiabatic regime of ^ = 2 scaling, 
the instanton-like optimal escape path can be found from 
linear equations of motion. 

We expect that the new <^ = 2 scaling and the scaling 
crossovers can be seen in various systems. Examples are 
modulated Josephson junctions, nanomagnets, and op- 
tically trapped Brownian particles, where escape in the 
presence of modulation has been already studied experi- 
mentally, albeit in different regimes. 

This research was supported by the NSF PHY-0071059 
and NSF DMR-0305746. 



APPENDIX A: VARIATIONAL EQUATIONS FOR 
THE ESCAPE PROBLEM 

We will consider here the optimal trajectory that the 
system is most likely to follow in escape, qopt(i)j ^'^d 
the random force fopt(i) that drives the system during 
this motion. The trajectories qopt(i), fopt(0 provide the 
absolute minimum to the functional 7?.[q, f] ©■ The vari- 
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ational equations have the form 



dt'J'{t-t')i{t') - X{t) =0, (Al) 



A= V(A K), 

q- K(q; - f(t) = 



(A2) 
(A3) 



(here and below we use the hat to indicate a matrix; 

v = a/aq). 

From Eqs. and HA1|) it follows that the optimal 
noise realization is expressed in terms of the Lagrange 
multiplier A and the matrix of the noise correlation func- 
tions 1^ (jSJ as 



^opt 



D- 



dt'^{t~t')\{t'). 



(A4) 



Therefore, for the optimal path, the activation energy 
functional © is 



1 1 



dtdt'\{t) ■ (fit - t')\{t') 



(A5) 



(note that the noise intensity D drops out from 7?., be- 
cause the correlation functions fij (t) are proportional to 
the noise intensity themselves). 

From the structure of the functional HA5|I (integration 
over time goes from — oo to cx) and the integral is non- 
negative) and from the fact that the system is initially 
in the vicinity of a stable state one immediately derives 
the boundary condition JSJ for t — > — oo. The arguments 
generalize to periodically driven systems the analysis of 
Ref. for stationary systems. 

Close to a periodic stable state q_a{t) time evolu- 
tion of X{t) can be described using the matrix /ia(t) = 
(dKi/dqj) ^, where the derivatives are evaluated for the 
state qa(i)- Because of the periodicity of qa(i), the ma- 
trix jlait) is also periodic in time. It determines the mon- 
odromy matrix 



Mait)=Ttexp 



t+TF 



dti ila{ti) 



where Tt is the operator of chronological ordering (cf. 
Appendix C). The matrix Ma shows how the distance 
between a point q and the cycle qa(i) varies over the 
modulation period in the absence of noise for small |q — 
qa(t)|: 

q(t + tf)- qa(t) = Ma{t)[q{t) - qa(t)] 

From the condition that c\a{t) is a stable state, the eigen- 
values of the matrix Ma are all less then 1 in absolute 
value: in this case the distance between q and qa(i) de- 
creases with increasing time. 

It is seen from Eq. ljA2|l that the monodromy matrix 
for A is the inverse transposed of Ma- Therefore its eigen- 
values are all larger than 1 in absolute value. Hence, if 



the system is in the stable state qa(i) for t ^ — oo, then 
\{t) ^ for t ^ -oo, and from ijMJ f{t) -f 0, too. 

For the periodic saddle qb(t) on the boundary of the 
attraction basin, one of the eigenvalues of the correspond- 
ing monodromy matrix exceeds 1 in absolute value. It is 
such a saddle-type boundary that can merge with an at- 
tractor at a saddle-node bifurcation we are interested in. 
Respectively, one eigenvalue of the matrix that describes 
time evolution of A is < 1. If A is pointing along the 
corresponding eigenvector, it will decay as t — *■ oo. Then 
f(i) will decay, too, from Eq. IjAip . This means that the 
system may asymptotically approach a saddle-type state. 
Note that there are no optimal paths that would go from 
one stable state to another, because the condition A ^ 
for i — > oo is not satisfied there. This explains the bound- 
ary condition (|SJ) for t — > oo. 

Because the function K is periodic in time, Eqs. (jAljl - 
(|A3|) with boundary conditions (|HJl have a periodic set 
of solutions. If q(t), f (t), A(t) is a solution, then q(t + 
TF),f{t+TF), X{t + Tp) is a solution, too. These solutions 
are heteroclinic orbits: they connect the states q_a{t), f = 
A = and qb(i),f = A = 0, which are also solutions 
of Eqs. ljAl|l - ljA3|l . Generally, only one heteroclinic orbit 
per period provides the minimum to the functional TZ lO . 



1. Escape in systems with a slow variable 

Eqs. ljAl|l - ljA3|l are largely simplified if one of the mo- 
tions in the system is slow and all other variables follow 
this motion adiabatically, i.e., their relaxation time i^"'' 
is much smaller than the relaxation time of the slow vari- 
able tr, at least for a part of the modulation period near 
a bifurcation point. We will assume that slow motion is 
described by qi ; the variable qi itself may be a periodic 
function modulated by a slowly varying factor Qi, as in 
the case discussed in Sec. Ill C. In this case we will be 
interested primarily in the factor Qi. 

Over a time ~ t^^\ the variables gi>i reach their equi- 
librium values q[{qi,t) for given Qi,t. They are deter- 
mined by the equations 



K,, 



,91 > 9, 



>!' 



(A6) 



In the absence of noise, ft — 0, the solutions of these 
equations are periodic for given Qi. As we will see below, 
the terms fi here are small and give small corrections. 

It is important that the periodic solutions q[yi with 
fi>i = be stable. From this condition and Eq. (|A2() it 
follows immediately that Ai>i = 0, otherwise the com- 
ponents Ai>i would exponentially grow in time, leading 
to the onset of a large force that would drive the system 
away from the state (|A6|) with given Qi. 

The nonzero component of the Lagrange multiplier Ai 
is determined by the component of the optimal force /i . 
The latter should overcome the restoring force on the 
slow variable Qi a-nd drive it from the stable to the unsta- 
ble state. But the slowness of Qi means that the restor- 
ing force is small, and therefore the force /i should be 
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small, too. It is smaller than the force that would be 
needed to overcome the restoring force for the fast vari- 
ables by at least a factor tf''^ /t-c- Therefore Ai is small 
as well. 

For given Ai(i) we can find all components /i>i from 
Eq. ijEIl- They are all cx 1/ir, and therefore to lowest 
order in tf^ jt^ they can be disregarded in the solution of 
Eq. l|X6)l for gi>i. 

The problem of escape is therefore reduced to a one- 
dimensional problem for the variable gi, the force /i, 
and the Lagrange multiplier Ai. In Eq. (|A3|I for gi, the 
functions gi>i should be replaced by g-^^ calculated for 
/, = 0. 

Further simplification occurs if the noise spectrum is 
smooth. The analysis here is different for the cases of slow 
and fast modulation, i.e. whether w^^r is small or large. 
The case LOpt^ 1 is discussed in Appendix C. Here 
and in Appendix B we consider the case LOpt^ ^ 1. The 
major effect of noise on the slow variable q\ comes from 
the noise spectral components at frequencies w ^ 1/ir. If 
the noise spectrum is fiat for such frequencies, the noise 
can be assumed to be white on the "slow" time scale. 
In other words, the correlation function </?ii(t) can be 
replaced by 2D5(t) [in this situation it is convenient to 
choose D from this condition rather than to define it by 
Eq. Q]. Then J^ii(t) = ^S{t). 

For a ID system driven by white noise of intensity D, 
Eqs. HAl|l - HA3p have a solution 

/i(t) = 2Ai(t)=gi-ifi, (A7) 

This reduces the variational problem of finding the opti- 
mal path to the known formulation for white-noise driven 
systems poj. 

APPENDIX B: REDUCED EQUATION OF 
MOTION FOR SLOW DRIVING 

In this Appendix we derive simplified equations of mo- 
tion for the case of slow driving where the relaxation 
time tr ^ Tp and the motion can be described in the 
adiabatic approximation. We will consider the vicinity 
of the adiabatic bifurcation point q — 0,t = 0, A = A^'^. 
A convenient basis for q is provided by the set of the 
right eigenvectors of the matrix fi = (dKi/dqj), where 
the derivatives are evaluated at the adiabatic bifurcation 
point. In this basis the equation of motion Q has the 
form 

qi ~ + ^ , K^.jkqjqk + K.-aSA'"'^ + K^.tt 

+ \K,,ut^ + Y K,.,,tqjt + Mt). (Bl) 

Here Ki.j^. = d^Ki/dqjdqk, Ki-A — dKi/dA, etc, with all 
derivatives calculated at the adiabatic bifurcation point. 



and SA'^'^ = A — A'^'^ . Since the function K depends on t 
only in terms of the modulation phase — ujpt, we have 
Ki;t oc ujp. The expansion in t in ljBl(l is, in fact, an 
expansion in cupt. 

Because the eigenvalue /ii is equal to zero, for small |q| 
relaxation of qi is much slower than relaxation of qi^i- 
In the absence of noise, over the relaxation time tr°^ the 
variables qiyi approach their quasistationary values for 
given qi and (f) — cupt. They can be obtained from the 
equations of motion ljBl|) for qi^i in which qi and fi{t) 
are disregarded, 

q, « -fi-' [K,,tt + K,,iiqj + K^.aSA-" + ...] (z > 1). 

(B2) 

Here, the major term is linear in copt. The full expres- 
sion is a series in qi,u!pt, and SA, and the omitted terms 
are of higher order in these variables. Because of the 
noise, qi^i will additionally perform small fluctuations 
with amplitude cx D^^^. 

The dynamics of the slow variable qi on times exceed- 
ing t^^^ is given by Eq. IjBip with i = 1, in which qiyi are 
replaced with their quasistationary values. To lowest or- 
der in qi, 6A^'^ , and ojpt only the linear in t term should 
be kept in Eq. HB2|) . This gives 

qi - aqf + (3 6A-^ ~ a^^ujpt)^ + h{t), (B3) 
a = ^Ki-ii, (3 = Ki-A, 

with 

Note that the coefficient 7 is independent of ti;^- 

Eq. (|B3|I reduces the multi-dimensional problem of 
random motion near a bifurcation point to a one- 
dimensional problem. In the case of a one-dimensional 
overdamped system driven by an additive periodic force 
it was derived earlier in Refs. Isil Is^ Besides the noise 
term, it differs from the equation of motion for stationary 
systems in the vicinity of a saddle-node bifurcation point 
[Tl| in that it has a term cx {uptY . 

Depending on the sign of SA^'^ , Eq. (jB3p has either two 
adiabatic solutions 

(gi ):% = Tsgn(a) [- {f3/a)SA^^ + {jcoptf] (B5) 

or none. For concreteness, we assume that the adiabatic 
solutions exist for SA^'^ < 0, i.e. af3 > 0. The solutions 
are even functions of time. They touch each other at i = 
for SA^'^ — 0. We assume that the periodic adiabatic 
states <i^'^i^{t) exist for all times provided SA'^'^ < 0. 

The term Ki-tt in Eq. (|Bip has to be equal to zero, 
otherwise the bifurcation point will be far from SA'^'^ = 
t = 0. On the other hand, the equation of motion 
(jB3|) may contain the term Cqiuspt, where C is a sum 
of Ki-it and appropriately weighted products Ki-nKi-f 
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This term can be eliminated by a linear transforma- 
tion 91 — >■ (?! + Cujpt/ 2a and renormalizations SA'^'^ — > 
(5^^'^ + CwF/2a/3, 72^72 + {C/2af. The renormalized 
7^ should be positive, if the stable and unstable adiabatic 
periodic states touch each other only for M^'* — t = 
and only once per period. 

The term cx qiujpt does not arise in the important 
case where the modulation is performed by an additive 
periodic force F{t), see Eq. Here, the adiabatic states 
'ia'^bi^) correspond to the minimum and maximum of the 
potential Uo{q) — F(i) • q, cf. Fig. ^ They merge first 
with increasing modulation amplitude A when the field 
component |_Fi| is at its maximum over t. This means 
that dtJi ~ 9q9tK = at the bifurcation point. 

As explained in Sec. Ill, the typical relaxation time 
near the bifurcation point does not exceed ti — 
{ajajp)^^^^ ■ If the correlation time of the noise fi{t) 
is much less than ti,t'^'^ and the power spectrum of the 
noise does not have singular features for high frequencies, 
then the dimensionless noise 

/» = (7c.j.)-ViW {T = t/ti) (B6) 

is effectively (5-correlated as a function of the "slow" time 
T, with (/(t)/(0)) = 2DSiT). From Eqs. ®, ^ the 
effective noise intensity is 

/oo 
dt^nit). (B7) 

APPENDIX C: REDUCED EQUATION OF 
MOTION FOR FAST DRIVING 

In this Section we consider the case where modula- 
tion near the bifurcation point is effectively fast, so that 
u!pt^ ^ 1. Here, throughout the modulation cycle the 
stable and unstable states c[a,b{t) stay close to each other 
and to the critical cycle qc(i) into which they merge at 
the bifurcation point A = Ac. Therefore the equation of 
motion can be expanded in 5q = q — qc(t), A — Ac, 
leading to Eq. (|22() . The expansion coefficients are peri- 
odic in time. 

It is convenient to start the analysis by simplifying the 
part of Eq. ^ 

6q = flSq, fi = fj.{t + Tp), 

that describes motion in the linear approximation in (5q. 
We introduce the matrix k{t,ti) such that 

k{tA) = Ttexp ^ dtifl{ti)j , (CI) 

where Tt is the operator of chronological ordering. This 
matrix satisfies the equation dk{t,ti) /dt = jl{t)k{t,ti) 
and gives the monodromy matrix M, 

M{t) = M{t + Tp) = k{t + Tp,t). 



The eigenvalues of the matrix M determine the 
evolution of 5q(t) in linear approximation. Over the pe- 
riod Tp, the coefficients of the expansion of Sq{t) in the 
right eigenvectors ey(t) of M change in times (we 
use Greek letters to enumerate eigenvalues and eigenvec- 
tors; they should be distinguished from the vector com- 
ponents, like Qi). The eigenvalues M^, are independent 
of time because of periodicity of M{t),fi{t). They are 
simply related to the Floquet exponents for the periodic 
state qc(i). 

At the saddle-node bifurcation, where stable and 
saddle-type states coalesce, one of the eigenvalues (for 
example. Mi) becomes equal to 1, whereas |Mi,>i| < 1. 
This means that the system is attracted to q_c{t) in all 
directions except for the critical direction ei(t); the dis- 
tance from qc(i) along ei{t) does not change over the pe- 
riod, in linear approximation. In what follows we choose 
ei to be real. 

For small SA of an appropriate sign, the state 
q_c{t) splits along ei{t) into a stable and an unsta- 
ble state. The system approaches the vicinity of these 
states along the directions e^^i over a short time 
Ti? max [l/|ln |Af^>i| |] ti^\ In contrast, the motion 
along ei is slow. 

The ei-component of Sq is the soft mode. We are 
interested in its dynamics on times long compared to 
t'f\Tp. The analysis is simplified by the fact that, for 
t — ti ^ [ri?/ In |My>i| I , the matrix k{t,ti) projects any 
vector Sq{ti) on the vector ei{t). In particular, 

k{t,U)ei{ti) « Kn{t,U)ei{t). (C2) 

This is a consequence of the transitive property k(<, ti) — 
k{t,t')k{t' ,ti) and the fact that, for an arbitrary vector 
(5q, we have M"(t)(5q Cei{t) for n 00. The func- 
tion Kii in Eq. (|C2|) is given by the expression 

Kii{t,ti) = ei{t) ■ k{t,ti)ei{ti). 

Here, ei is the left eigenvector of the matrix M, which 
corresponds to the eigenvalue Mi = 1, and we use nor- 
malization §1 (t) ■ei{t) = 1. The matrix element kh (i, ti) 
is periodic, Kii{t + Tp, ti) = Kii{t, ti). 

Equation of motion for Kii(t, ti) for large t — ti follows 
from Eq. 

d 

— [Kii(t,ti)ei(i)] = ^iii{t)Kii{t,ti)ei{t), 

fin{t) = ei{t)- fi{t)ei{t). (C3) 

Close to the bifurcation point, the component of (5q 
along the vector ei (t) has a slowly varying factor. In con- 
trast, the components of Sq along the vectors e^^i are 
"fast" . Over time ^ ti^^ they reach quasiperiodic values 
for a given value of the slow component, and then fluctu- 
ate with amplitude cx D^^^. From (|22|l . the quasiperiodic 
values are quadratic in the slow component and therefore 
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small. As a consequence, the slow motion is indeed one- 
dimensional, 

<5q(t)«Qi(i)«:ii(i,ii)ei(0. (C4) 

The instant ij here is arbitrary; Qi{t) contains a multi- 
plicative factor that depends on ti (but Sq{t) is indepen- 
dent oiti). The time ti drops out of all final expressions, 
see Sec. IV. 

The equation for Qi(t) is obtained by substituting 
Eq. ((nil) into Eq. ^ and then muhiplying Eq. ^ 
by the vector ei{t) from the left. This gives 

=/C(Qi,i)+ei(i)-f(i), (C5) 
IC{Qi,t) = ^Kl,{t,ti)Ql{e,{t) ■ vfK^ + <M(a^i^i), 

where Ki ei • K. 

In the absence of noise, the solution of Eq. HC5|I is a 
sum of smooth and oscillating parts, Qi{t) — Q^"^{t) + 
Qosc^ The term Q^™ remains nearly constant on the 
time scale tf, whereas Q°^'^ ^ ujpQ°^'^. It is seen from 
Eq. (jnUl that Q°^'= cx 6A. The term Q™ is much 
larger. An equation for Q**™ can be obtained by aver- 
aging Eq. (|C5|) over time. It has the form 

Q™i = a'(Q™)2+/3'M + /'(t). (C6) 

The coefficients a' , (3' in Eq. ljC6|l are given by the 
expressions 

a' - i(Kn(t,ii)(ei(t)- V)'i^i)^, 

/?' = {K-l{t,t;)dK,ldA)^, (C7) 

where ^-^^ means period-average centered at time t, 

{g)^^Tp^ r^''\t' g{t\ti). (C8) 

Jt-TF/2 

The result of the averaging ljC8|l is independent of t for 
time-periodic as in the case of the coefficients a', (3', 
and therefore a',/3' are independent of t. 



The function f'{t) in Eq. (|C6|I is a random force, 

f'{t)^K^,\t,U)e,{t)-i{t). (C9) 

Eq. I|('6|l has the same form as the equation for the 
soft mode in the adiabatic limit (|12|l in the absence of 
the term c>c (cupt)'^. For a'P'SA < the system has a 
stable and an unstable stationary solution Q™. These 
solutions are given by an equation similar to Eq. H19(l . 

QT, = Tsgn(a')(-/3' 5A/a')'^' (CIO) 

(in what follows without lost of generality we set a' > 0). 

Typical values of Q^"^ are oc |(i4p/^, as seen from 
Eq. ljC10|l . They largely exceed the amplitude of the fast 
variables in 6q, which are all oc &4, in the neglect of noise. 
The relaxation time of Q''™ is U = |2aQ™|-i oc \SA\-^^^. 
It is much larger than tf close to the bifurcation point. 
The condition LUpti- 3> 1 was the major approximation 
made in the derivation of Eqs. IjCSp . I|C6|) . besides the 
condition of the weak noise. 

A transformation from Q''™, t to reduced variables 

Q = a'i/2Q--, T = a'^^h (Cll) 

allows us to write Eq. (jC6p in the compact form (|23|) . The 
random force /(t) = f'{a'^^^t) is effectively (5-correlated. 
From Eqs. ©, HC9I) . its intensity is 

/oo 
-OO 

xei(i + h) ■ ifih - t2)ei{t + h))^. (C12) 

Here, if> is the matrix of the noise correlation functions 
®. 

As a result of the period-averaging over i, in Eq. I|C12|) 
the integrand becomes a function of ti ~ t2, and there- 
fore the integral over ti is independent of t2. Still, it 
depends on ti, but this dependence will drop out of the 
final expressions for observable quantities, in particular 
the activation energy of escape. 
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